Entanglement evolution of two qubits under noisy environments 
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The entanglement evolution of two qubits under local, single- and two- sided noisy channels is 
investigated. It is found that for all pure initial states, the entanglement under a one-sided noisy 
channel is completely determined by the maximal trace distance which is the main element to 
construct the measure of non-Markovianity. For the two-sided noisy channel case, when the qubits 
are initially prepared in a general class of states, no matter pure or mixed, the entanglement can 
be expressed as the products of initial entanglement and the channels' action on the ma^dmally 
entangled state. 
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I. INTRODUCTION 

In realistic quantum information processing [l|, 
quantum systems are open to essentially uncontrol- 
lable environments that act as sources of decoherence 
and dissipation^, Therefore the study of entangle- 
ment evolution of quantum systems taking into account 
the effect of environment has become more and more 
interesting [5l-[To|. In Ref.Q Konrad et al. provided a 
direct relationship between the initial and final entangle- 
ment of two qubits with one qubit subject to incoherent 
dynamics. It has been shown that given any one-sided 
quantum channel, the concurrence of output state corre- 
sponding to any initial pure input state of interest can 
always be equivalently obtained by the product of the 
concurrence of input state and that of the output state 
with the maximally entangled state as an input state. 
However, for the two-sided quantum channel case or the 
mixed initial state case, the product of the two concur- 
rences only provides an upper bound for the concurrence 
of interest. It is obviously important to find the exact re- 
lations of that for two-sided noisy channel case and that 
for the mixed initial state case. 

Over the past decades, the conventionally employed 
Markovian approximation with the assumption of an in- 
finitely short correlation time of environment has expe- 
rienced more and more challenges due to the advance of 
experimental techniques 3] . Recent studies [ll| - |25j have 
shown that non-Markovian quantum processes play an 
increasingly important role in many fields of physics. In 
order to quantitatively study the non-Markovian dynam- 
ics, following the first computable measure of Markovian- 
ity for quantum channels introduced in Ref.[2()|, some 
measures for the degree of non-Markovianity have been 
introduced [13, HI]. It is very interesting to investigate the 
explicit relationship between the entanglement evolution 



and the non-Markovianity. 

In this paper, we consider the entanglement evolution 
of two qubits under local noisy channels. The main aim of 
this work is to analyze if and to what extent can the fac- 
torization law given in Ref.[^] be generalized. As a cen- 
tral result, we demonstrate that for the two-sided noisy 
channels, even when the qubits are initially prepared in 
some mixed states, the concurrence can be expressed as 
the products of the initial concurrence and the channels' 
action on the maximally entangled state. 

The paper is organized as follows: In Sec|lll we present 
the model and its analytical solution. In Sec lIIII we give 
the entanglement dynamics of two qubits under local, 
single- and two- sided non-Markovian channels. In Sec lIVI 
two applications are given. Finally we give a conclusion 
of our results in Sec El 



II. THE MODEL 

We consider a system formed by two nonintcracting 
parts A and B, each part consisting of a qubit s ^ a,b 
locally interacting respectively with a reservoir Rs = 
RatRb- Each qubit and the corresponding reservoir are 
initially considered independent. The dynamics of each 
part, consisting a qubit with excited state |e) and ground 
state \g) which is coupled to a reservoir of field modes 
initially in the vacuum state, can be represented by the 
reduced density matrix Q 



'Plmhs{t)\' 



(1) 



in the qubit basis {|e), \g)}, where the superscript S — 
A,B represents part S. The function hs(t) {h{t) for 
short) is defined as the solution of the integrodifferen- 
tial equation 



f{t-h)h{t), 



(2) 



'Electronic address: 
t Electronic address: 



jungl@bit.edu.cn] 



zoujian@bit.edu.cn] 



with f{t — ti) denotes the two-point reservoir correlation 
function which can be written as the Fourier Transform 



2 



of the spectral density J(a;): 

f{t — ti)~ / dujJ{uj)exp[i{LOo~u!){t~ti)]. (3) 



The exact form of h(t) thus depends on the particular 
choice of the spectral density of the reservoir. It should 
be noted that the dynamical map of the model can be 
Markovian (iff \h{t)\ is a monotonically decreasing func- 
tion of time) a.nd non-Markovian (iff \h(t)\ increases at 
some times) [2^. 

Next, we are ready to use the reduced density matrix 
elements given in Eq. ([T]) to construct the reduced density 
matrix for the two-qubit system. Following the procedure 
presented in Ref. 25], we find that in the standard prod- 
uct basis B = {jl) = |ee),|2) = |e5),|3) = l^e), |4) - 
\gg), }, the diagonal elements of the reduced density ma- 
trix p(t) for the two-qubits system can be written as 



Piiit) = piimhA{tr\hB{t)\^ 

P22{t) = \hA{mP22{o) + pii(o)(i - ihBimi 

pz,{t) = \hB{mpzm + pii(o)(i - \hA{t)?% 

pn{t) = 1 - (pn(t) + p22{t) + P33{t)), 

and the nondiagonal elements are 

Pl2{t) = Pl2mhA{t)\''hB{t), 
Pisit) = Pl3{0)hA{t)\hB{t)\'', 
Pu{t) = pu{0)hA{t)hB{t), 
P23{t) = P23{0)hA(t)hB{t), 

P2i{t) = hA{t)[p2m + pi3(0)(l - \hB{t)\^)]. 
P3i{t) = hB{t)[p^m + P12(0)(1 - \hA{t)\'')l 



(4) 



(5) 



and Pij{t) = Pjiit)- We should note that if we let 
hA{t) — h{t) and hB{t) — 1, this reduced density ma- 
trix represents the dynamics of the two qubits in the 
single-sided channel case. 



III. ENTANGLEMENT EVOLUTION 

In this section, we give the entanglement evolution for 
single-sided channel case {hA{t) = h{t), hB{t) = 1) and 
two-sided channels case for pure and mixed initial states. 
Firstly, we assume that the qubits are initially prepared 
in the pure state 



I*) 



ci|ee) -I- C2\eg) + Csl^e) -|- c^lgg), 



(6) 



where Ci{i — 1,2,3,4) are complex and satisfy the nor- 
malization condition X]i=i I'^iP ~ 1- The elements of 
the density matrix corresponding to this initial state can 
be written as Pij{0) = CiC*,{i,j = 1,2,3,4). In what 
follows, we consider the single-sided channel case. Sub- 
stituting hA{t) = h{t), hsit) = 1 and Pij{0) = CiC* into 
Eqs. (|H) and ([5]) we obtain the diagonal elements 



Pii(i) 

P22{t) 
P33{t) 
Piiit) 



clci\h{t)\^ 

C*2C2\hit)\^ 

c5c3 + ctci(l-|Mi)P) 

clC4 + 4C2{1 ~ \h{t)\'') 



(7) 



and the nondiagonal elements 

Pl2(i) =Cic5|/l(t)|2, 

Pi3{t) = ciclh{t), 
Pii{t) = cicXh{t), 

P23{t) = C2Clh{t), 
P2i{t) = C2CXh(t), 
P3i{t) = C3C| -I- CiClil 



(8) 



and p,j{t) = p*,{t). 

To quantify the entanglement we use Wootters 
concurrence [30|. From Eqs. ([7]) and we find that 
the concurrence of the two qubits is 



C*(t) = C*(0)|/i(t)|, 



(9) 



where Ciii(O) = 2|ciC4 — C2C3I is the initial entanglement 
of state © . Equation ^ shows that the entanglement 
can be factorized as two terms whose physical meaning is 
given as follows. We consider the extreme case in which 
the two qubits are initially prepared in the maximally en- 
tangled state, i.e., C3-(0) = 1. Then equation ^ reduces 
to Cii>(t) = |^(i)|, that means that \h{t)\ stands for the 
time evolution of the entanglement for the maximally en- 
tangled initial state under the single-sided channel. Then 
we can say that the entanglement reduction under single- 
sided noisy channel is equal to the product of the initial 
entanglement C>ii(0) and the time evolution of the entan- 
glement for the maximally entangled initial state. That 
is, the dynamics of the entanglement is completely deter- 
mined by the time evolution of the entanglement for the 
initial state with maximal entanglement. 

The factorization law given in Eq. ^ is valid for single- 
sided channel case. One may ask to what extent it can 
be generalized to the two-sided channels case. In what 
follows, we will show that this can be done by constrain- 
ing the initial state to a certain class of states. We also 
let the two qubits be initially prepared in the pure state 



®. Substituting py(0) 



into Eqs. (g]) and © and 



after some calculations we find that the concurrence at 
time t is given by 



C*(i) = max{0,g(t)}. 
The explicit expression of Q{t) is 

Q{t)^C^{0)\hAmhB{t)\X, 



(10) 



(11) 



where X = 1 ~ |cip/|c2C3 — ciC4\y^ with £, = {I — 
|/iA(i)P)(l - l^s(i)P)- Equation shows that for the 
case of two-sided noisy channels, the entanglement evo- 
lution is not equal to the product of the initial entangle- 
ment C^(0) and the maximal entanglement evolution of 
single-sided channel |ft.A(i)| and \ hB{t)\ but an additional 
term X should be taken into account. We should note 
that if for some time intervals A" < 0, the entanglement 
sudden death occurs. Interestingly if we consider the case 
ci = 0, that is, the two qubits are initially prepared in 
the state 



1$) 



C2\eg) +C3\ge) +C4\gg), 



(12) 
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with x;,=2 |c» 

reduces to 



1, we find that X = 1 and Eq.fTOI 



C^{t) ^ C^{0)\hA{t)\\hB{t)\, 



(13) 



with C$(0) = 2|c2C3| being the entanglement of initial 
state (HH). Equation (|13l) indicates that for the pure ini- 
tial state p^ . the entanglement evolution under two- 
sided noisy channels is equal to the product of the ini- 
tial entanglement C$(0), (the maximal entangle- 
ment evolution under single-sided channel A) and 
(the maximal entanglement evolution under single-sided 
channel B). 

As we know, quantum mechanics and quantum infor- 
mation processing are not constrained to pure states. So 
we give the generalization of the factorization law to the 
mixed state. We consider the case in which the qubits 
are initially prepared in the mixed state 



P(0) 



/ 

b z e 

z* c f 

V e* /* d . 



(14) 



where 6, c and d are real numbers and satisfy 6 -I- c -I- 
d = 1. For this state it is easy to check that the initial 
concurrence Cp(0) = max{0, \/bc+ \z\ — \Vbc— \z\\} and 
the terms Pij{0) in Eqs.(|l]) and ([S]) can be written as 

pii(0) = 0, pi2(0) = 0, pi3(0) = 0, pu{0) = 0, 

P2l(0)-0, P22(0)-6, P23{0) = Z, P24(0)=e, 
P3l(0) = 0, p32{0) = Z*, P33(0) = C, P34(0) = /, 
P4l(0) = 0, p42(0) = e*, P43(0) = /* P44(0) = d. 

After some calculations we find that the concurrence at 
time t is given by 



Cp{t) ^ Cp{Q)\hAmhB{t)\. 



(15) 



That is the factorization law for the two-sided noisy chan- 
nel case. We should note that the initial states ([12]) and 
p4p have a common property, that is, the number of ex- 
citations in the system is not more than one. We denote 
these states NOE states. Then we can conclude that 
for the two-sided noisy channels case, when the qubits 
are initially prepared in NOE states, no matter pure or 
mixed, the concurrence can be expressed as the prod- 
ucts of initial concurrence and the channels' action on 
the maximally entangled state. 



IV. APPLICATIONS 

In this section we give two applications of our results. 
The first one is the entanglement preservation. Equa- 
tions ([T(7)) and show the direction for entanglement 
preservation, that is, choosing appropriate parameters to 
make the values of |^A(i)| and \hB{t)\ around 1. As an 
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FIG. 1: Dynamics of entanglement as functions of "fot. 



example, we investigate the detuning case of a Lorentian 
spectral density [J 



1 



27r [ujQ 



A)2 + A2 



(16) 



where A = ujq — uJc is the detuning of ujc and wq, and uJc is 
the center frequency of the cavity. A defines the spectral 
width of the reservoir and is connected to the reservoir 
correlation time tji = X~^. 70 is related to the decay of 
the excited state of the qubit in the Markovian limit of 
a fiat spectrum, the relaxation time scale T5 over which 
the state of the system changes is then related to 70 by 
T5 = 7(7^. We evaluate the reservoir correlation function 
f(t — ti) using the spectral density J{uj) and obtain 

/(i-ii) = i7oAexp[-(A-zA)(f-ii)]. (17) 

Solving equation ([2]) with this correlation function, we 
find 



h{t) 



-i(A-iA)t 



coshi — 
2 



A - zA . , fdt 

; — smh — 

d \ 2 



(18) 

with d = Y^(A - iA)2 - 270A. From Eq.dTH]) we find that 
the value of \h{t) \ is determined by the parameters A and 
A. For the two-sided channel, we denote A^ and Xa for 
channel A and A^ and Ag for channel B. 

Figure [1] shows the dynamics of the entanglement as 
functions of 70^ for different cases. In FiglU Ci is the en- 
tanglement evolution of the initial state ip = l/^/2{\eg) + 
\ge)) for the single-sided channel case and C2 is that for 
the two-sided channel case. We choose the parameters 
Aa = Ab = 0.570 and Xa — Xb — O.OI70 to make 
the values of |/ia(OI and |ft.B(t)| around 1. Clearly, the 
entanglement is protected well. 

In FiglU we also give the entanglement evolutions for 
two-sided channel case of initial state ip = l/^/2{\eg) + 
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\ge}) (C3) and initial state (p = l/^/2(|ee) + \gg)) (C4). 
Here, we choose the parameters Aa = As = and 
= Xb = O.OI70. Comparing C3 and C4 we find that 
for initial state ip (NOE sate), the entanglement is al- 
ways large than zero except some discrete time instants. 
This can be understood as follows. For the NOE sates 
the factorization law ()13p is valid. The entanglement is 
determined by and |/is(i)| of which the values are 

not less than zero for all the parameters. However, for 
the initial state (j) (not NOE sate) , the entanglement sud- 
den death and the entanglement sudden birth can occur. 
This can be understood from Eqs.([TU| and PT|) that in 
this case X can be less than zero for some time intervals. 

The second application is to show the connection of the 
entanglement and the non-Markovianity of the channels. 
The measure for non-Markovianity defined in Ref. (27j is 

7V= max / dta{t), (19) 

Pl,2(0) J^yo 

where a{t) is the rate of change of the trace distance 
which can be defined as 

a{t) = j^D{p,{t),p2it)), (20) 

with D(pi, P2) = P2I being the trace distance of the 

quantum states pi and p2- Here, D{pi,p2) describes the 
distinguishability between the two states and satisfying 
Q < D < 1. We should note that if there exists a pair of 
initial states and a certain time t such that a[t) > 0, the 
process is non-Markovian. Physically, this means that 
for non-Markovian dynamics the distinguishability of the 
pair of states increases at certain times. This can be in- 
terpreted as a fiow of information from the environment 
back to the system which enhances the possibility of dis- 
tinguishing the two states. Following Ref. [27], reference 
(311 1 presented a practical idea for directly measuring the 
non-Markovian character of a single qubit coupled to a 
zero-temperature bosonic reservoir. It is easy to check 
that any pair of initial states satisfying the conditions in 
the theorem given in Ref. ^31.] definitely owns the same 
maximal trace distance 

D{p,{t),p2{t))^\h{t)l (21) 

namely, \h{t)\ stands for the maximal trace distance of 
one part [A or B) of our model. From Eq.([9]) we can 



find that the entanglement evolution under single-sided 
noisy channel is equal to the product of the initial en- 
tanglment C^(0) and the maximal trace distance |/i(i)|. 
We also know that \h{t)\ also stands for the channels' ac- 
tion on the maximally entangled state. It is evident from 
this relation that any increase of the entanglement im- 
plies occurring of the non-Markovian dynamics. Hence a 
conceptually simple way to quantify the degree of non- 
Markovianity of an unknown quantum evolution would 
be to compute the amount of entanglement at different 
times within a selected interval and check for strict mono- 
tonic decrease of the quantum correlations. This mea- 
surement has been given in the very recent paper [2^. 

V. CONCLUSION 

In summary, we have studied the entanglement dynam- 
ics of two qubits under local, single- and two-sided noisy 
channels. We find that in the case where the initial state 
is pure, and only one of the subsystems undergoes a non- 
Markovian channel, the equation of motion for entangle- 
ment presents the form of a simple factorization law -the 
second term contains only information about how the 
maximal entanglement is affected by the dynamics, and 
the first term scales the second by the initial amount of 
entanglement. For more realistic scenarios, where both 
parts are influenced by the local environments, when the 
qubits are initially prepared in NOE states, no matter 
pure or mixed, the concurrence can be expressed as the 
products of initial concurrence and the channels' action 
on the maximally entangled state. Using these factoriza- 
tion laws we have found the way to protect entanglement 
from decay and given a connection of the entanglement 
and the non-Markovianity of the channels. 
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